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TEMPERATURES, THERMAL STRESS, AND SHOCK IN HEAT- 
GENERATING PLATES OF CONSTANT CONDUCTIVITY AND 
OF CONDUCTIVITY THAT VARIES LINEARLY WITH 
TEMPERATURE 
By S. V. Manson 
SUMMARY 

Working formulas are presented for the steady-state temperatures 
and thermal stress in heat -generating infinite plates of constant con- 
ductivity, and of conductivity that decreases linearly with t emp erature 
as the temperature Increases , for the case in which the heat is generated 
uniformly throughout the plate thickness and Both faces of the plate are 
equally cooled. In addition to exact formulas for the variable- 
conductivity plates , simpler approximate formulas are derived, and cri- 
teria for their applicability are indicated. 

It is shown that, of all planes in the plate, the plate surfaces 
are always under the greatest tension, and the midplane is under the 
greatest concession. 

A criterion Is indicated for determining the surface cooling condi- 
tions under which the thermal shocks at the surface and midplane will he 
smaller than, equal to, or greater than the steady-state thermal stresses 
at those planes. The shocks exceed the steady-Btate stresses only when 
the initial surface-cooling conditions of the transient state are more 
severe than the surface cooling conditions of the steady state. The 
criterion is expressed for heat-generating plates by a simple relation 
between a steady-state heat-transfer parameter of the heat-generating 
plate and the Biot number of the transient state. 

The dimensionless parameters governing the transient temperatures 
and thermal stresses in materials of linearly varying conductivity are 
derived by a similarity study of the conduction equation and boundary 
conditions of the transient state. A numerical technique for solving 
the transient-state equations is indicated in detail. The method is 
employed to obtain numerical values of the transient temperatures and 
stresses for a variety of parameter combinations, selected to test the 
thermal-shock criterion. The results of the calculations substantiate 
the shock criterion. 
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As a result of the study an insight is obtained into the steady and 
transient temperature and stress mechanisms in heat -generating plates of 
constant conductivity and of conductivity that varies linearly with tem- 
perature. 


INTRODUCTION 

Thermal stresses in the steady and transient states have acquired 
increasing interest for aircraft: in ductile-material applications, as 

in turbine disks and blades; and in brittle -material applications, as in 
nozzle diaphragms, combustion-chamber linings, high- temperature coatings, 
and recently in solids that generate heat inter nally . 

A substantial literature exists on the steady-state thermal stresses 
arising from temperature nonuniformity in solids of cylindrical and rec- 
tangular geometries. The greater part of the literature on transient- 
state stresses (shock), however, is concerned with suddenly heated or 
cooled bodies of initially uniform temperature . No systematic study of 
shock in bodies of initially nonuniform temperature nor any criterion 
for ready prediction of when the shock will exceed the steady-state 
stress was found. 

An additional consideration which appears not to have been studied 
is the question whether thermal conductivity variations produce signifi- 
cant stress effects, in those applications for which the operating tem- 
perature range is large. Examples of solids of substantially varying 
conductivity are low-carbon steel, which halves its conductivity in the 
range from room temperature to 1450° F, and Inconel X which doubles its 
conductivity in the same temperature range (ref . l) . Cermets and other 
high-temperature materials may also exhibit substantial conductivity 
variations with temperature. If conductivity variations are considered, 
it is desirable to identify the parameters governing the temperature in 
both the steady and transient states, inasmuch as such an identification 
facilitates design and the choice of test conditions. 

A specific problem of current interest is that of thermal stress 
and shock in solids that generate heat internally. If the power per unit 
volume of solid is high, and if the solid has low conductivity or 
unfavorable creep properties, the thermal stress or shock may be exces- 
sive. For heat-generating solids the effect of substantially decreasing 
conductivity with increasing temperature is also of considerable interest, 
inasmuch as such conductivity variation aggravates the steady stresses. 

As part of its research in materials for aircraft, the NACA Lewis 
laboratory has investigated the steady and transient temperatures , and 
the steady and transient thermal stresses, in infinite plates which gen- 
erate heat uniformly throughout- their interior in the steady state and 
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which., on suddenly ceasing to generate heat, are simultaneously im mersed 
in an enviro nme nt of constant temperatur e . The study has been made for 
cases in which the plate conductivity is constant and for cases in which 
the conductivity decreases linearly with temperature as the temperature 
increases . 

The steady-state temperature distribution in a heat -generating 
plate of linearly varying conductivity was not found in the literature, 
and hence is derived herein; for comparison and completeness, the well- 
known steady-state temperature distribution in a heat-generating plate 
of constant conductivity is also presented. The transient temperature 
distribution in plates having an initially nonuniform temperature dis- 
tribution and a thermal conductivity that varies continually with loca- 
tion and time was also not found in the literature. Such a plate can- 
not- be readily analyzed by the existing procedures for composite plates, 
the conductivities of which are generally treated as uniform and as 
time-independent in a region of any one material; hence an extension of 
the well-known finite-difference procedures was employed. The transient 
temperatures in a constant-conductivity plate having an initially non- 
uniform temperature distribution are presented in an infinite-series 
solution in reference 2. For consistency with the calculations on 
variable-conductivity plates, and also to eliminate series-convergence 
considerations for the short time-periods of interest in the thermal- 
shock study, the finite-difference procedure was preferred over the 
infinite-series solution for the constant-conductivity plates. 

In the study of heat -generating plates reported herein, formulas 
for the steady-state temperatures and stresses are presented. A cri- 
terion is indicated for determining whether the steady stress or thermal 
shock will be the greater under anticipated conditions of operation. 

The dimensionless groups governing the transient temperatures and ther- 
mal stresses in materials of linearly varying conductivity are derived. 
A numerical technique is indicated for solving the transient-state con- 
duction equation for variable-conductivity materials, and typical curves 
of transient temperature and stress are presented. 


AKALTSIS 

General Stress Considerations 

The relation for thermal stress in a plate of thickness 2L 
(fig. l), in which the temperature varies only in the x-direction, 
and in which end effects are not considered, is given in reference 3, 
and in the notation of this report is 
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(Symbols are defined in appendix A- ) The geometric variable x may he 
changed to the d im ensionless form x/L, and the temperature variable T 
may be changed to 8 defined by the relation 


8 = T - T e 


( 2 ) 


where T e is a constant. If (3E is satisfactorily constant over the 
temperature range of interest in an application, and if use is made of 
the facts that 



equation (l) may be written 



From equation (2), the 8 value at the plate surface is 


®s “ Tg - Te 


( 4 ) 


and if the steady-state surface tengperature is denoted by 


T 


s,o> 


B.O 


T 


S ,0 


- T c 


( 5 ) 
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and & B>0 is a constant. 

If a is employed to represent both cry and a z , which are equal 
to each other, and if o* is defined as 



equation (3) may be rewritten in the dimensionless form 

o*= a -to.). = + i f 1 d (l 

^0s,o \®s,o/ 2 / \0s,o/ \i 




Equation (7) involves no explicit time dependence; if 0/0s,o is known 
throughout the plate the dimensionless stress is determined regardless 
of whether 0/0 8j)O is a steady or transient value. The first term of 

the right member of equation (7) is the local value of e / d B)Q at a 

plane chosen for consideration, and may vary as x/L varies; the second 
term is the average value of 0/0 s ,o over the plate thickness, and for 
an existing temperature distribution is constant . The third term' is a 
"moment of a temperature-moment" and may vary with x/L, but if the tem- 
perature distribution is symmetrical about the plate center line, 

(x/L = 0), the term vanishes, and the local dimensionless stress is equal 
to the difference between the plate average and the local values of 
0/08,0* Thus in the case of a symmetrical temperature distribution the 
thermal stress arises simply from the tendency of the local plane to 
ejfpand or contract to a length different from that corresponding to the 
plate average temperature. 

The case of a temperature distribution symmetrical about x/L = 0 
is the case under consideration in this report. As has been mentioned, 
equation (7) reduces, in this case, to 


a* = . g to) . = -LU) + toj) 

f3E0 s ,o \0s,o/x/L \0s,o/av 


( 8 ) 



6 


MCA TO 2988 


Equation (8) indicates that there is no stress at the plane x/L where 
the temperature is the same as the plate average temperature. If a 
positive value of o* denotes tension, and a negative value of cr* 
denotes compression, equation (8) indicates that planes at which the tem- 
perature exceeds the plate average temperature are in compression, and 
planes at which the temperature is less than the plate average tempera- 
ture are in tension. Thus, for any plate temperature distribution sym- 
metrical about the midplane the surface and midplane are under the 
greatest stresses of opposite sign. For heat -generating plates equally 
cooled at both surfaces, the surface is under the greatest tension, and 
the midplane under the greatest compression, of al 1 planes in the plate. 


Steady-State Temperature and Stress in Heat-Generating Plates 

The actual magnitudes of surface and midplane stresses of heat- 
exchanging plates can be determined if the distribution of Q[9 q , 0 is 

known explicitly. For heat-generating plates in steady state, the fol- 
lowing temperature distributions are derived in appendix B: 

(a) For a plate material of constant thermal conductivity. 


G 


s,o 


» 1 + 


q m Ii 
2K0 


s.o 


t - (if] 


(9) 


(b) For a plate material having a conductivity that varies linearly 
in the temperature range between environment and surface temperatures. 



(An approximate formula for ®/0s,o i s discussed in appendix B.) 

Steady-state stresses for heat -generating plate of constant con- 
ductivity. - Hie average value of 9]Q s , q with equation (9) is 
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or 




/ q'" l2\ 

{^oj 


(1 


On substituting equations (9) and (ll) into equation (7) or (8), 



Equation (12) indicates zero stress at — = — — , where the temperature 

h V3 

is equal to the plate average temperature . At the plate surface, 

» ± 1 , the stress is 


a* _ °x,o( 1 - v ) 2 / g"' L 2 \ 

' - eE9 S;0 3 \2K0 S )0 J 


at the plate mi dp lane, x/L a 0, the stress is 


a * _ a m,o(l" v ) 1 ( q 1 " L 2 \ 
m, ° = W Q s,o 83 3 \2K» S>0 J 


(i: 


( 1 < 


Steady-state stress for heat -generating plate of linearly varying 
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and if only positive values of A are considered (positive values of 
A correspond to plates of conductivity that decreases as temperature 
increases , as indicated in appendix 33), 



By use of equations (10) and (16) in equation (8), the stress at any 
x/L can he computed. At the surface, x/L = -tl, the stress is 



at the midplane, x/L = 0, the stress is 



(Approximate formulas for a*^ Q and are derived in appendix C.) 

Transient Temperatures and Thermal Shock 

The 0/Qs,o distribution (and hence the stress) in the transient 
state is governed hy the familiar equation 


po K “ 3J ( K S 


(19) 
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in. tlae plate interior; and at the surface, "by the standard relation 

*" s (s). = (20) 

Equation (19) contains no source term because it is assumed for the 
analysis herein that the hbat -generating plate has suddenly ceased to 
generate heat at the instant that it is immersed in the transient cooling 
environment. It is noted that, whereas in the steady state the constant 
T e is arbitrarily assignable, • the use of 8 S in equation (20) indicates 

that in the transient state the value assigned to T e must be the con- 
stant tenperature of the environment into which the plate is immersed on 
ceasing to generate heat. In practice it is convenient to take T e the 
same for both the steady and unsteady states, even though the steady- 
state coolant may not have the temperature T e . 

A calculation procedure for obtaining numerical values of 0/® S}O 

from equations (19) and (20) is indicated in appendix D, and in appen- 
dix E a derivation is presented of the dimensionless groups governing 
Q/9 B) o in the transient state. In terms of the derived groups, values 
of 0 /Q s Q computed for an arbitrarily prescribed initial temperature 
distribution may be presented in a form useful for similar systems. 

From a knowledge of the 9/6 s ^ 0 distribution the transient thermal 

stress (shock) is obtained by use of equation (8) of this section, as is 
indicated In some additional detail in the Calculation Procedure of 
appendix D. 


CRITERION FOR DOMINANCE OF STEADY-STATE THERMAL 
STRESS OR TRANSIENT THERMAL SHOCK 

In establishing a criterion for whether the steady-state stress or 
the thermal. shock will be greater, the heat-transfer process in the 
steady state and in the early stages of the transient state may be exam- 
ined in some detail, and inferences be drawn concerning the simultaneous 
stress variation. The ultimate guide is equation (8) of the Analysis 
section, which indicates that the stress is proportional to (0 av - 9 x /i) • 

Inasmuch as the greatest stresses occur at x/L = 0 and ±1, that Is, at 
the plate midplane and surfaces, the following considerations are directed 
at a determination of the manner in which the midplane and surface tem- 
peratures vary In relation to the plate average temperature. 

For convenience in discussion the half -plate may be assumed divided 
into N thin segments of equal width 5. If 1 square foot of frontal 
area is considered, and the plate is assumed to generate q’" Btu per 
second per cubic foot uniformly throughout its interior In the steady 
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state, then each segment of width 6 generates q n ' 5 Btu per second. 
The steady-state temperature profile (i.e., local slope), is such that 
each segment receives all the heat generated hy segments more interior 
than itself and transmits not only the heat received hut also the addi- 
tional amount q 5 that it generates) hence the profile steepens out- 
ward. The segment hounded hy the plate surface receives (iT-l) 6q IM , 
generates 8q' M , and transmits NSq'" to the sink. If h s t represents 
the steady-state heat-transfer coefficient. 


hst (Ts ~ T c ) — hg-fc. (®s “ ®c) = = q'" L (2l) 


Associated with the steady- state temperature profile and heat flow is a 
thermal-stress distribution, everywhere proportional to the difference 
between the average temperature and the local temperature. 


s,o/av 


s,o 


where the subscript o denotes steady state. 


At the instant that the plate ceases to generate heat and is 
immersed into the coolant of temperature T e , the temperature profile in 
the plate interior is that of the steady state and the heat flow through- 
out the plate and into the plate surface is the same as in the steady 
state. Accordingly, in a short time interval At of the transient 
state each element of thickness 6 (except possibly the element next 
to the plate surface) delivers outward all the heat received from seg- 
ments more interior than itself, and in addition gives up 5q M, AT 
heat units from its own volume. Inasmuch as the heat source has ceased 
to generate, the amount of heat 8q'” At given up by each segment is 
derived from the enthalpy of the segment itself) hence the temperature 
of every interior segment falls by q " ' At/ pc degrees. The level of 
the interior temperature is thus decreased, but the temperature profile 
everywhere except near the surface is exactly as before. Accordingly, 
in the next time interval At each interior segment, delivering heat 
at the rate determined by the local temperature profile (slope), again 
gives up 5q ' " At heat units. It is seen that in a finite period of 
the transient .state, before the plate interior can be affected by the 
cooling conditions at the surface, the interior temperature level must 
fall uniformly with time, at the rate q" ’/pc at every point) and the 
interior temperature profile must remain for a period exactly the same 
as in the steady state, regardless of the nature of the surface cooling 
conditions. It is accordingly seen that the only influence the midplane 
region can exert on the plate average temperature early in the transient 
state is to promote a decrease of the average temperature at the rate 
q'" /pc. 
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The major factor capable of modifying the rate at -which the average 
temperature changes is the surface temperature, which is at all times 
responsive to the external cooling conditions. Hence the stress-shock 
criterion must emerge from a study of the temperature variations in the 
plate as determined by the three possible kinds of cooling condition at 
the plate surface, as follows: 

If h is the constant heat -transfer coefficient in the transient 
state, Ts^o the initial temperature of the plate surface, and T e the 
constant temperature of the transient-state environment, the heat flow 
from plate to environment at the instant of immersion is 


b(T S j 0 - ^e) = h0 s ^ o 
and three possibilities exist: 

( a ) k®s,o = h s .j. ( 9 Sf0 ® c ) 

O’) M s , 0 < i s t ( 8 s,o - 6 c) (23) 

(=> M S)0 > h, t (e s>0 - e c ) 

Case (a) . - In case (a) of equations (23) the plate surface is sub- 
ject to initial transient-state cooling conditions equal to those of the 
steady state. Then the environment removes heat at exactly the rate of 
its delivery from plate midplane to surface, with no heat retention or 
extraordinary heat flow occurring anywhere in the plate. Then not only 
all interior temperatures, but the surface temperature as well, initially 
fall at the rate q ,T, /pc, and for a period the temperature profile 
throughout the plate (except for second-order changes at the surface) 
will remain the same as in the steady state. During this period of con- 
stant temperature profile the average plate temperature of course also 
falls at the rate q M, /pc; hence for this transient-state period the 
thermal stresses everywhere in the plate must, to first-order accuracy, 
be constant at the values of the steady state. Ultimately the surface 
temperature decreases enough that h0 s is perceptibly smaller than 
h®s,o* At this time 


( 22 ) 


(ET - l) 5q m < h0 g < h0 s ^ o = NSq m 
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and 



(24) 


The decrease of the temperature slope at the plate surface indicated by 
inequality (24) arises because the segment next to the surface can no 
longer deliver heat to the environment at the rate prevailing in the 
steady state and must retain a fraction of the enthalpy it formerly deliv- 
ered from its own volume. Thus the temperature of the surface segment 
can no longer fall at the rate q'"/pc by which the temperatures of 
nearby points are falling. Inasmuch as the average plate temperature 
reflects the change at the surface, the average temperature now also falls 
at a rate slower than q " '/pc, although not as slowly as the temperature 
of the surface segment itself. At this time both the midplane and sur- 
face temperatures approach more closely than formerly to the average tem- 
perature, and the thermal stresses at the two planes must decrease to 
below the steady-state values. 

It is thus seen that if the cooling conditions at the surface at 
the beginning of the transient state are equal to those of the steady 
state, the stresses everywhere in the plate must for a period remain 
effectively constant at the steady-state values; and that the first 
changes in surface and midplane stresses are decreases to values below 
those of the steady state. 

In case (a) it is noteworthy, as for the surface segment just con- 
sidered, that every segment of the formerly-heat-generating plate is at 
some time subject to an undiminishing heat inflow at one face and a 
diminishing sink strength at the other face, and that the segment adjusts 
the rate of delivery of its own enthalpy to satisfy the heat-flow con- 
ditions at both faces. In this manner the plate temperature profile 
flattens inwardly -until the midplane is affected by the cooling conditions 
at the surface. Throughout the period of initial profile flattening the 
stresses at the surface and midplane decrease. 

As the surface temperature decreases, the ability of the environment 
to remove heat also decreases; and inasmuch as the steadily decreasing 
sink strength of the environment always requires time to affect the 
average and midplane temperatures, the temperatures of the interior will 
always drop more rapidly than the surface temperature. Thus in case (a) 
the surface and mddplane stresses will decrease steadily to zero. 

Case (b). - In case (b) the plate initially delivers heat to the 
surface at a rate greater than the environment can remove; hence the tem- 
perature of one or more segments near the surface rises immediately on 
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immersion, while simultaneously the interior temperatures fall at the 
rate q ’ " /pc . The average temperature must fall at a rate slightly 
lower than q"’/pc; hence the midplane and surface temperatures "both 
come closer to the average than in the steady state, and "both stresses 
must immediately decrease. For the same reason as in case (a), that the 
average and midplane temperatures each feel and respond to the changing 
sink strength of the environment later than the surface does, the tem- 
perature profile must flatten steadily inward, and the midplane and sur- 
face stresses must decrease steadily to zero. 

Case (c). - In case (c) of relation (23) the environment removes 
more heat than the plate interior initially delivers to the surface. 

The excess heat derives from segments near the surface, and the surface 
temperature must fall at a rate exceeding q’"/pc. Hence the average 
plate temperature will also fall more rapidly than q m /pc, hut less 
rapidly than the surface temperature. The midplane temperature, as 
previously discussed, falls initially at q'"/pc. Hence both the surface 
and midplane stresses must immediately rise to values greater than those 
of the steady state. 

Inasmuch as the effect of the surface cooling conditions must pene- 
trate into a substantial portion of the interior before the average tem- 
perature can begin to fall at the same rate as the surface temperature, 
the surface stress will continue to rise for a period. Similarly, the 
rate of decrease of midplane temperature lags behind that of the average 
temperature, and the midplane stress will also continue to rise. Ulti- 
mately the diminishing sink strength of the environment as surface tem- 
perature decreases, coupled with increasing heat flow from the plate 
interior, will lead to a leveling off and subsequent steady decrease of 
the surface stress. At a later time the midplane stress will reach its 
peak and subsequently subside. 

Thus, it is seen that for formerly-heat-generating plates it is 
only in case (c), h0 s ^. > h g ^. ( 9 8j0 - ® c ), that the surface temperature 

initially falls more rapidly than the interior temperatures; and only 
in this case do the shocks at the surface and midplane exceed the steady- 
state stresses. (It is of interest that even in case (c) there are 
intermediate segments for which the stress initially decreases, because 
the stress becomes zero whereever the local temperature equals the aver- 
age; in case (c) the location of the average temperature moves closer 
to the surface on plate immersion.) 

Linearly varying conductivity. - The foregoing discussion has been 
for plates of constant conductivity. The effect of linear increase in 
conductivity as temperature decreases is now briefly considered. In 
general, the local temperature profile may be expected to dominate over 
the conductivity variation in determining the heat flow. Hence the 
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general “behavior of variable-conductivity plates may be expected to be 
similar to that of constant-conductivity plates. If it is noted, however, 
that AK/K per degree decrease in temperature is greater in the interior 
than near the surface, it follows that in cases (a) and (b) the tempera- 
ture profile will flatten more rapidly, and the midplane and surface 
stresses will decrease more rapidly, than if conductivity is constant. 

In the case (c) the surface stress will always tend initially to 
increase to a value above that of the steady state. The initial stress 
variation at the midplane, however, depends on the degree to which 
h0 s ^ o exceeds h s t (0s, o - 0c)> because the favorable conductivity var- 
iation at the midplane may permit the temperature there to decrease at 
a rate equal to or exceeding that of the plate average temperature . 
Inasmuch as the surface stress always tends to increase when h0 s ^ o 

exceeds h B ^. (0 S Q - 0 C ), a calculation of the transient stresses will 
be necessary. For a fixed ratio h0 s ^ o /h s ^. (0 g ^ o - 0 C ) greater than 

unity, the maximum surface and midplane shocks in the variable- 
conductivity plate will be a smaller multiple of the steady-state values 
than in the constant-conductivity plate; the favorable K variation in 
the plate interior reduces the rates at which both the surface and mid- 
plane temperatures depart from the plate average temperature. 

Shock criterion in terms of hl/he and q ,M L 2 /2Ke0 8 , o - ~ The con- 
dition of equation (23a), 


h0 s,o = *st ( 0 s,o - 0 c) ( 23a ) 

which is the condition separating dominance by steady-state thermal 
stress from dominance by transient thermal shock, may be expressed as 
a relation between the transient-state parameter hl/Kg and the steady- 
state parameter q' n L 2 /2Ke0B,o of heat-generating plates, as follows: 
From equation (2l) 

( s 8,0 - S C> - 1"' L 

Then when equation (23a) is satisfied. 


h0 s , o = 


( 25 ) 


On multiplying both sides of equation (25) by L/Kg0 s ^ o 
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hi _ 2 ( Z ' Z JL.) 

Ke " \Z£e e s,o) 


( 26 ) 


Thus if equation (26) is satisfied the steady-state stress and maximum 
transient shock will he equal to each other. The shock will exceed the 
stress only if 


f >-2 

K e 



(27) 


RESUITS AMD DISCUSSION 
Steady State 

Temperatures . - Formulas for the steady-state temperature distribu- 
tions in heat -generating plates of constant conductivity, and of con- 
ductivity that changes linearly with temperature, are derived in appen- 
dix B and have been indicated in the Analysis section (eqs. (9) and (10)). 
The temperature variable employed in the derived formulas is 0, defined 
as T - T e , where T is the local plate temperature, and T e is a 
constant reference temperature. (T e is not necessarily the temperature 
of the steady-state cooling environment. The use of 0 puts the steady- 
state temperatures into the same form as the one useful in describing 
the plate transient temperatures that would occur on sudden exposure of 
the plate to a cooling environment of temperature T e . ) In the derived 
formulas the temperatures are expressed in dimensionless form by use of 
the ratio 0/0 3 , Q , a parameter appearing among the dimensionless groups 

shown in appendix D to govern the transient temperatures in plates of 
linearly varying conductivity. From the ratio Q[Q q ^ q the quantity' 

[jp/ 0 s,o)“ l] is readily obtained and, being equal to (T - T gjO )/0 g ^ o , 

gives a direct measure of the temperature difference between the surface 
and an interior point in the heat -generating plate.. 

The formulas derived in appendix B show that in heat -generating 
plates of constant conductivity the ratio Q[Q e}0 depends on x/L and 

on q ,M L 2 /2K0 8 0 ; in heat -generating plates of linearly varying conduc- 
tivity the ratio 0/0 g ^ o depends on x/L, q ,M L 2 /2K e 0 g ^ o , and Kg^/Kg. 
Curves illustrative of the distribution of ®/®s,o 8X6 presented in 

figure 2 in the form {(©/0 S , o) " 3 against x/L. Parts (a) to (d) of 
figure 2 are for the four values 0.0025, 0.005, 0,010, and 0.020 of the 
parameter q"' L^/2Kg0 SjO , and in each part of the figure plates of con- 
stant conductivity, and plates of K g ^ 0 /Kg values of 0.5 and 0.2, are 
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represented, (it is noteworthy that low values of K s ^ 0 /K e do not- 

necessarily imply a large slope of conductivity with temperature . A 
moderate conductivity slope with temperature, combined either with a 
large temperature difference & B}0 ) or with a low value of reference con- 
ductivity Kg, or both, can result in values of Kg^ 0 /K e substantially 

lower than unity. ) The curves of figure 2 reflect the relations derived 
in appendix B: 


For plates of constant conductivity the quantity is 

parabolic with x/L, and at any x/L increases in direct proportion to 
q ,M L 2 /2K0 S ^ O . For plates of linearly varying conductivity, and for 

appropriate combinations of q'"L 2 /2K e 0 s Q and Kg Q /K e as discussed 
in detail in appendix B, the quantity |(0/0s,q) - lj is also essentially 
parabolic with x/L, and at any x/L increases very nearly directly with 
q ,n l£/2Kg0 6 ^ o and inversely with Kg^ 0 /Kgj these characteristics are 

illustrated in figure 2(a) to (c) wherein all the plates are behaving 
very nearly as constant-conductivity plates, with conductivity equal to 
that at the .temperature of the plate, surface. Figure 2(d) illustrates. 


however, in the case 



0.2, that for combinations of q ,M L 2 /2K e 0 8 ^ o 


and Kg^ 0 /k e outside the suitable range, as discussed in appendix B, 

variable- conductivity plates differ to a significant degree from constant- 
conductivity plates. A measure of the accuracy with which the constant- 
conductivity relation approximates the variable-conductivity relation is 
discussed in appendix B, where it is also shown that if K decreases as 
temperature increases there is a limiting value beyond which the param- 
eter q 1 " L^/2K e 0 g ^ o cannot increase. 


Stresses. - Relations have been indicated in the Analysis section 
(eqs. (13), (l4), (17), and (18)) and in appendix C for the steady-state 
stresses in heat -generating plates of constant conductivity, and of 
conductivity that decreases linearly with temperature as temperature 
increases. Note has also been made in the Analysis section that in heat- 
generating plates equally cooled on both surfaces the greatest tensile 
stress occurs at the plate surfaces and the greatest compressive stress 
occurs at the plate midplane . Values of dimensionless steady-state 
surface stress a*^ Q , and mLdplane stress cr*^ 0 , in heat -generating plates 

of constant and of linearly varying thermal conductivity are shown in 
figure 3 for values of q ,M L 2 /2K 0 0 S ^ O from 0 to 0.0145 and for values 

of K s ^ 0 /Kg from 0 to 0.15. Figure 3 reflects the stress formulas 
derived In the Analysis section and in appendix C: 
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For plates of constant conductivity "both. o B)0 and <% 0 increase 
in direct proportion to q' M L 2 /2K0 S ^ Q . For plates of conductivity that 

decreases linearly with temperature the stresses are at first very nearly 
directly proportional to q m L 2 /2K e 0 S and inversely proportional to 

K s ^o/K e - For the combinations of q" ’ L 2 /2K e 0 s ^ o and Kg^/Kg for 

which these direct and inverse proportionalities occur the plate is 
behaving as one having a constant conductivity equal to that at the tem- 
perature of the plate surface. The range of q n ' L 2 /2K e 0 s ^ o in which 

cr* is linear with q' 11 L 2 /2Kg0 B ^ o depends on the value of K s ^ 0 /K e , 
and the range increases as K g ‘ increases. Thus for K g Q /K e of 

0.15, curvature away from a straight line begins to occur at 

q’" L 2 / 2K e 8 g Q of 0.003, while for Kg^/K^ above 0.375 linearity 

exists over the range of q ,M L 2 /2Kg0 s ^ o (0 to 0.0145) shown in figure 3. 
For each Kg^/Kg smaller than unity there is a value of q"' L 2 /2K e 0 s ^ o 
above which the variation of a * ceases to be linear; whenever such cur- 
vature occurs the inverse proportionality between a* and K B ^ 0 /K e 

also ceases to be valid, (it is noteworthy that when cr* iB linear with 
qin cr* is directly proportional to plate thickness, provided 

that h s t( 0 s,o “ ®c^ the heat dissipation per unit surface area, is 
constant. The constancy of h ai .(0 g ^ o - 0 C ) implies that, as plate 

thickness increased, the internal heat -generation rate q ,M is decreasing 
as l^Ii. ) 


Transient State 

Governing dimensionless parameters. - If a plate of arbitrary 
steady-state temperature distribution is suddenly immersed in a cooling 
environment of constant temperature T e the subsequent temperature and 
thermal stress histories of the plate depend on the plate initial tem- 
perature distribution; on T e ; on the transient-state heat-transfer 

coefficient; on plate thickness; and on plate thermal diffusivity, which 
may vary with temperature. In order that a systematic study and presen- 
tation of transient temperatures and stresses may be made, a knowledge 
of the parameters that govern the transient state is necessary. 

General transient-state parameters for variable-conductivity plates: 
In appendix E it is shown that if the conductivity of a plate material 
varies linearly with temperature in the range between plate and environ- 
ment temperatures, the parameters governing the transient -temperature 
ratio 0/0 s ^ o at x/L (and hence the dimensionless stress cr* at 

x/L), are hl/K e , a e T /l 2 , and K s ^ 0 /K e . The first two parameters are 

the Biot and Fourier numbers, respectively, with plate thermal conduc- 
tivity evaluated at the constant temperature T e of the cooling envi- 
ronment. The parameter Kg 0 /K e is the ratio of plate conductivity at 
the initial surface temperature T s ^ 0 to the plate conductivity at 
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temperature T 0 . The indicated three parameters, then, are suitable for 
the systematic study and presentation of transient temperatures in plates 
of linearly varying conductivity. (For a fixed value of K 3 , 0 /K q the 

Biot and Fourier numbers can equally well he based on Kg^ 0 instead of 
KeJ hence hL/Kg^ 0 , a g ^ q t /L^ , and Kg Q /K e could equally well be regarded 
as the governing parameters . ) 


It is noteworthy that although the thermal conductivity of a plate 
material 1 may vary substantially with temperature the constant conductiv- 
ities K^ or Kg^ 0 in the Biot and Fourier numbers, together with the 

constant ratio K g 0 /K e , are completely adequate to define the transient 
temperatures, for a prescribed initial temperature distribution. It is 
also noteworthy that plate materials having different conductivities at 
T e , and different conductivity slopes with temperature, will have equal 


transient values of 0/0 g , o at x/L if the values of hL/K e , a e T/L 2 , 
and * 8 , 0 / K e are equal in the two plate systems (provided, of course. 


that the initial distributions of dfQ B ^ Q with x/L are also equal. 

The requirement of similar steady-state temperature distributions sig- 
nifies, for example, that plates having uniform initial 0/ 0 B ,o can 

compared only with other plates having uniform initial 0/0 g ^ o ; while 
plates having linear, or parabolic, initial 0/0 g ^ o distributions can 

only be compared, respectively, with other plates having linear, or 
parabolic, initial distributions). 


Parameters for heat -generating plates: The derived steady-state 

formulas show that for heat -generating plates of constant conductivity 
the steady-state value of 0/0 g ^ o at x/L depends on q M ’ L 2 /2K0 g ^ o , 

and for heat -generating plates of linearly varying conductivity the 
steady-state value of 0/0 s , o at x/L depends on q , "L 2 /2K e 0 s ^ o anh 
K g 0 /k g. Thus if a plate generates heat uniformly in the steady State, 

and then ceases to generate heat and is simultaneously immersed in a 
cooling environment of temperature T e , the transient-temperature ratio 
®/ d s,o a ^ X A- ds completely determined if a specification is made of 
q’" L 2 /2K e 0 g ^ o , Kg 0 /K e , hl/K e , and . (For plates of constant 

conductivity, q"' L^/2K e 0 g ^ o is, of course, the same as q ,M L 2 /2K0 0 ^ O , 
and K g ^ 0 /K e equals unity.) 


Criterion for Dominance of Thermal Stress or Shock and Typical 

Variations of Transient Temperature and Stress 

In the Analysis section it was indicated that the relative magni- 
tudes of steady-state stress and transient shock depend on the relative 
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magnitudes of h0 SjO and h s t (0s,o - ®c)j or equivalently, for heat- 
generating plates, on the comparative magnitudes of hL/Kg and 
q ,n L^/2K e 0 s ^ o . The results of temperature and stress calculations for 

several values of the stated parameters are presented in figures 4 to 9; 
figures 4 to 6 are for plates of constant conductivity, and figures 7 to 
9 are for plates of Kg^ 0 /K e equal to 0.151. In al 1 cases the plates 

are assumed to have generated heat uniformly throughout their interiors 
in the steady state, and are assumed to have ceased suddenly to generate 
heat and to have been immersed in a cooling environment of constant tem- 
perature T e - All plates are assumed to he cooled equally on both faces, 

by constant heat-transfer coefficient h' between plate and cooling 
environment . 

Constant- conductivity plates, hL/K = 2(q ,n L^/2K0 SjO ). - Figure 4 

is for plates of constant conductivity, for a steady-state value of 
q'" I?/2KQ b , q equal to 0.01, and for a transient-state value of hL/Kg 

equal to 0.02. In figure 4(a), 9/9 3}0 is shown against x/L for several 

values of ax/l? . For the calculation an interval of slightly less than 
0.02 was employed for arcfl?, as satisfying the convergence criterion 
indicated in the Calculation Procedure of appendix D; most values of 
<xt/l2 are omitted from the figure. In figure 4(b) the dimensionless 
surface stress 0 * and the dimensionless midplane stress 0 * axe shown 

against a x/L 2 . 

The temperatures and stresses in figure 4 illustrate the conclu- 
sions obtained in the Analysis for plates of constant conductivity, for 
hL/K = 2 (q' M L^/2K0 S ^ O ): The temperature profile, and hence the stresses 

throughout the plate, remain for a period effectively constant at the 
steady-state values; and after the period of constancy the stresses at 
surface and midplane decrease. 

Constant-conductivity plates, hL/K < 2(q'" l2/2K0 s ,o) • - Figure 5 
is for plates of constant conductivity, for a steady-state value of 
q" ' J?/2K Qq }0 equal to 0.01 and for a transient-state value of hL/K 

also equal to 0.01. Figure 5 illustrates the conclusions of the Analysis 
for the constant-conductivity case, hL/K < 2(q"' l2/2K0 s q ): The tem- 

peratures near the surface initially rise, while the interior tempera- 
tures fall uniformly; hence the surface and midplane stresses immediately 
decreaae to values below those of the steady state. The figure also 
shows that the temperature profile throughout the plate ultimately adjusts 
itself so that every segment transmits not only al 1 the heat it receives, 
but also gives up very nearly the same amount of its own enthalpy as 
given by the other segments. When this adjustment occurs the temperature 
throughout the plate falls very nearly uniformly with time, as indicated 
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"by the curves for axf'S? from 0.40 to 0.48. Inasmuch as this condition 
is lik e the one previously presented in figure 4 (hL/K = 2q ,M L 2 /2K0 g ^ o ), 

the t emp erature and stress behaviors subsequent to <xt/L 2 of 0.48 will 
closely resemble the behavior of the previously presented case. 

Constant-conductivity plates, hL/K > 2(q' " L 2 /2K©s,o) • - Figure 6 is 

for plates of constant conductivity, for a steady- state value of 
qitt j2/2K0 s q equal to 0.01, and for a transient-state value of hL/K 

equal to 0.03. 

Figure 6(a) shows that initially the temperatures near the surface 
fall more sharply than those of the interior, reflecting the fact that 
the coo ling conditions at the surface are more severe than in the steady 
state. The temperature curve steepens progressively inward until, for 
the case shown, ax/j? equals about 0.38, after which the shape of the 
temperature curve is such that each segment gives up very nearly the same 
enthalpy as every other segment; thereafter the temperature falls very 
nearly uniformly throughout the plate, until ax[l? equals at least 
0.937. The subsequent temperature history will closely resemble that 
for the case Mj/k = 2(q ,M L 2 /2K0 S ^ O ), illustrated in figure 4. 

Figure 6(b) shows that at the instant of plate immersion o* and 
both increase, of more rapidly than a*. At cct/l 2 of about 0.4, 
for the case shown, both stresses have slowed their rate of increase 
almost to zero. The stress a* is effectively at its peak value, 0.0099, 

a 50 percent increase above the 0.0066 value of the steady state. At 
oct/L 2 of 0.937, 0* has passed its peak value and has started to decrease, 

while o* is still increasing very slowly prior to attaining its own 
peak, which will differ only slightly from the value at anr/L 2 of 0.937. 
Thus the maximum transient value of 0 *, like 0 *, is about 50 percent 

higher than the steady-state value. 

The temperature and stress variations shown in figure 6 illustrate 
the conclusion in the Analysis that when hL/K exceeds 2(q ,u L^/2K9 8 f0 ) 
the surface temperature will initially fall more sharply than the uni- 
formly falling interior temperatures, and hence that the surface and mid- 
plane stresses will increase initially to values above those of the 
steady state before decreasing to zero. 

Variable-conductivity plates, hL/K e = 2(q ,M L 2 /2K e 9 B}Q ). - Figure 7 

is for plates of conductivity that decreases linearly as temperature 
increases; the conductivity variation is represented by ttie ratio 
K s ,o/&e = 0.151. The figure corresponds to a steady-state value of 
qiu x J 2 /2K e 0g 0 equal to 0.01 and a transient-state value of hl/K equal 
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to 0.02. In figure 7(a), 9/d 3}Q is shown against' x/l for several 
values of a^r /l? . For the calculation an interval of 0.0756 was used 
for aer/l 2 , as required by the convergence criterion indicated in the 
Calculation Procedure of appendix D. The solid curves represent tempera- 
tures computed with the transient-state conductivity variation taken into 
account. The dashed curves represent temperatures computed when, with 
the actual steady-state temperature distribution of the variable- 
conductivity plate taken as the starting distribution, the conductivity 
variation in the transient state was neglected and the initial surface 
conductivity was employed for the entire transient- state calculation. 
Figure 7(b) shows dimensionless surface and midplane stresses, o* and 

a m , respectively, and the solid and dashed curves have the same signifi- 
cance as indicated for figure 7(a). 

Figure 7(a) shows that the steady-state temperature distribution 
(ctetr/L 2 = 0) is much steeper than the steady-state distribution shown 
in figure 4(a) for constant-conductivity plates. The difference in 
steepness for equal values of q ,M L 2 /2K e 0 SjO arises from the much lower 

plate surface and interior conductivities when K Q}Q /& e equals 0.151 
than when K s ^ 0 /K e equals unity . Figure 7(a) shows that although the 

surface cooling conditions at the start of the transient state are equal 
to those of the steady state, the interior temperatures fall more rapidly 
than the temperatures near the surface. The difference in cooling rates 
of the interior and surface regions arises because although the conduc- 
tivity increases everywhere in the plate as the temperature falls, the 
value of AK/K per degree temperature drop is greater in the interior 
than near the surface. 

The dashed lines of figure 7(a) show that the computation based on 
surface conductivity, a value not much in error for the plate region near 
the surface, but noticeably higher than the interior conductivities, 
leads to an initial overestimation of the heat flow from the plate 
interior, and to a corresponding overestimation of the early rate of 
temperature drop in the plate interior. From the relation h0 8 = q" 

it is seen that when the net heat flow is overestimated 9 e must be 
higher than when q" is accurately evaluated. Hence, for the time span 
under consideration (a e T/L 2 S 1.2), the surface temperatures on the 

dashed curves are higher than those of the more accurate solid curves. 

The comparative long-term temperature histories, after the conductivity 
everywhere in the plate exceeds the initial surface value, would require 
further calculation to establish. 

Figure 7(b) shows that as in the case of constant-conductivity 
plates when hL/K = 2(q m L 2 /2K0 S ^ O ), the surface- stress is for a short 

period effectively constant at the steady-state value, and then slowly 
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decreases. (it may "be noted that the steady-state stress is much higher 

than for the constant-conductivity plate at the same value of 

q' ’* J?f?K Q Q Q}0 , because of the steep temperature variation in the case 

K s ,o/ K e = 0-151, as previously indicated. The steady-state stress is 
slightly lower than the value shown in figure 3 for K s } 0 /K e = 0.150 

because of the slightly better conductivity represented by 

K Sj0 /Ke = 0.151.) Uni ike the case of constant-conductivity plates, the 

midplane stress decreases immediately upon plate immersion, reflecting 
the approach of midplane and average temperatures as the midplane con- 
ductivity and temperature drop increase. Tor values of OgT/Ij greater 
than 1.2, both a* and a* will slowly but steadily decrease to zero. 

The dashed curves of figure 7(b) show that because of the over- 
estimated flattening of the temperature curve exhibited in figure 7(a), 
the stresses at both surface and midplane are initially underestimated 
if computations are based on constant conductivity equal to that at the 
initial surface temperature. Thus for this case, at a e r f 1? of 1.0, 

the surface stress is underestimated about 8 percent, and the midplane 
stress about 20 percent. Neither underestimation is very significant, 
however, if the steady-state stresses are tolerable, because as indicated 
in the Analysis, and as indicated also by figure 7(b), the transient 
stresses at surface and midplane do not exceed the steady-state values 
when hL/Ke equals 2(q ,M L 2 /2K0 Sj Q ) . (A conservative estimate of the 
stresses would result if the less laborious computation with constant 
conductivity were based on conductivity equal to that at the initial 
midplane temperature . ) 

Variable-conductivity plates, hi JKq < 2(q ,M L 2 /2K e 0 s Q ). - Figure 8 
is for plates of conductivity variation represented by K B ^ Q /K e » 0.151. 
The figure corresponds to values of q”' L 2 /2K e 0 s ^ o and hl/lCg both 

equal to 0.01. As a whole, the transient temperatures (fig. 8(a)) -and 
the transient stresses (fig. 8(b)) exhibit behaviors very similar to 
those in the analogous case of constant-conductivity plates indicated in 
figure 5. The modifications introduced by the conductivity variation 
are basically the same as those already discussed in figure 7. Also, 
as in figure 7, the use of constant conductivity equal to the initial 
surface value, instead of the variable local values, initially over- 
estimates the cooling rate of the plate interior and underestimates the 
cooling rate of the surface region? hence underestimates the early sur- 
face and midplane stresses. 

Variable-conductivity plates, hh/Ke > 2(q''' L 2 /2Ke0s,o)- - Figure 9 

is for plates having a linear conductivity variation represented by 
K s ^ 0 /k e = 0.151. The figure is for q ,M L 2 /2K e 0 s ^ o equal to 0.01 and 
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for hL/K e equal to 0.03}' thus hL/Kg > 2(q m L 2 /2K e G B Q ) and. the 

initial transient-state surface cooling conditions are more severe than 
those of the steady state. For the temperature and stress calculations 
shown in figure 9, an interval of slightly less than 0.0747 was used for 
a^-r/L . The solid and dashed curves have the same significance as indi- 
cated in figure 7. 


The solid lines of figure 9(a) show that in the early transient 
state (cL e x/Tu smaller than about 0.3) the midplane and surface regions 
cool more rapidly than the intermediate plate region. The midplane 
region cools relatively rapidly, because AK/K per degree temperature 
decrease is relatively high. The surface region cools rapidly because 
the environment extracts heat at a rate greater than can be supplied 
by the plate interior, the net heat flow from which is primarily governed 
by the temperature profile and only secondarily by the conductivity var- 
iation. By the time a e T / L 2 is about 0.3, however, the plate tempera- 
ture profile has adjusted itself so that each segment cools at very 
nearly the same rate as every other segment, a transient-state adjustment 
which was also noted in the cases of constant -conductivity plates herein 
considered. The dashed curves of figure 9(a) show, as in previous figures, 
that the initial surface conductivity overestimates the net heat flow 
from the plate, mainly from overestimating the cooling rate of the mid- 
plane region. 


Figure 9(b) shows, as in the case of constant-conductivity plates, 
that when hL/K > 2(q I " L /2K e 0 s ^ o ) the surface stress rises on immer- 
sion of the plate into the transient-state environment. The surface 
stress reaches a maximum and then decreases. It is noteworthy that the 
maximum percentage rise of surface stress for K s ,o/&e of 0-151 is 
about 20 percent of the steady-state value, compared with the 50-percent 
increase noted for constant-conductivity plates having the same 
q n, L^/2K e 0 s ^ o and hl/K e (fig. 6). The favorable effect of an increas- 
ing conductivity on the transient-state stresses was indicated in the 
Analysis. As in previous figures, the constant-conductivity calculation, 
represented by the dashed curves of figure 9(b), underestimates the 
surface stress, but for i^he conditions of figure 9 the effect is less 
than 5 percent of the exact value. 


The midplane stress is seen in figure 9(b) to undergo an initial 
decrease, because of the excess local-temperature decrease as compared 
with the average plate-temperature decrease. The midplane stress sub- 
sequently increases, however, to a value about 10 percent above that of 
the steady-state value, as compared with approximately 50-percent 
increase noted for 4 in the constant-conductivity plates of figure 6. 
(At Oex/L 2 of about 1.0 in figure 10, the midplane stress is still 
Increasing very slowly, but its achievement of a maximum and subsequent 
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decrease are immi nent.) The constant-conductivity calculation underr 
estimates the midplane stress hy about 18 percent of the exact value and 
changes the time at which the stress starts to rise after its initial 
decrease. As previously indicated, advantage nay he taken of the smaller 
labor involved in the constant-conductivity calculation, while obtaining 
a conservatively high approximation of the stresses, by employing the 
conductivity corresponding to the initial temperature of the plate mid- 
plane . 

Typical variations. - It is seen from figures 4 to 9 for plates of 
constant conductivity, and of conductivity which decreases linearly as 
temperature increases, that the three possible variations of early 
transient-state temperature and stress are accounted for by exploration 
of the combinations (a) hl/Ke =j 2 (q* 11 l 2 /2K e 0 SjO ), 

(b) hL/K e < 2(q'" l 2 /2K e 9 s ^ 0 ), and (c) hl/Kg > 2(q" ' l 2 /2Kee s f 0 ) . Hence 
it is seen that the comparison of hL/Ke and 2(q” 1 lP/ZK e 8 BjQ ) is a 
suitable criterion for determining whether the steady stress or the 
transient shock will be the greater. It may also be seen that whereas 
values of q* " l^/2Ke0s,o and hl/Ke other than those herein considered 
will lead to magnitudes of ®/® s ,o ahd cr* different from those shown 

in the figures, the nature of the transient-state variations will be 
similar to those of either case (a), (b), or (c). That is, the curves 
of figures 4 to 9 are typical cooling and shock curves for former ly- 
heat -generating plates of constant conductivity, and of conductivity that 
decreases linearly as temperature increases. 

Conclusions . - From the discussion of figures 2 to 9 the following 
conclusions may be drawn: 

1. The steady-state temperatures and stresses increase as 
q 1 " L^/2Kg0 SjO increases and as K B}0 /Kq decreases. 

2. In the early transient state the deep interior of the plate cools 
at a rate determined entirely by the steady-state local temperature pro- 
file and by the local conductivity variations as the early cooling pro- 
ceeds; that is, the cooling rate at the midplane early in the transient 
state is independent of the cooling conditions at the surface. The sur- 
face region, on the other hand, cools at a rate determined simultaneously 
by the heat flow from the interior, the strength of the external sink, and 
the local conductivity variation. The parameters governing the transient 
temperature ratio S/0s,o at any point x/L are q ,,, L 2 /2K e 0 B „, K s ^ 0 /K e , 

hl/Ke, and Oqt/iP. 

3. The severity of the thermal shock, as compared with the steady- 
state thermal stress, may be satisfactorily evaluated by a comparison 
between the initial surface cooling condition of the transient state 
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h0 g ^ o and the steady-state surface cooling condition h st( S s,o ~ e c)* 

For formerly -heat -generating plates as herein described, the stress- 
shock compar ison may he made, equivalently, by a comparison of the Biot 
number hL/k e and the steady-state quantity 2(q ,M . The 

thermal shocks at both midplane and surface exceed the steady-state stress 
only when h0 g ^ o exceeds b gt (0 s ^ o - 0 C ), or equivalently, when hi/Kg 

exceeds 2(q'" L2/2Kg0 SjO ). 

4. If the conductivity of the plate material increases linearly as 
the temperature decreases, the maximum shocks of plate surface and mid- 
plane, when hL/K e > 2‘(q"' L 2 /2Kg0 g ^ o ), are a smaller multiple of the 

steady-state values than when the conductivity of the plate material is ' 
constant . 

5. Use of the initial surface conductivity for calculation of the 
transient temperatures and stresses initially overestimates the cooling 
rate and can substantially underestimate the early transient stresses. 

A conservative constant-conductivity calculation would be obtained by 
use of the steady-state conductivity of the plate midplane. 

6. The conclusions derived from figures 4 to 9 for the transient 
state substantiate the indications obtained from the analysis presented 
herein. 


RESULTS AND CONCLUSIONS 

1. Working formulas are derived for the steady-state temperatures 
and thermal stresses in heat -generating infinite plates of constant con- 
ductivity, and of conductivity that decreases linearly with temperature 
as the temperature increases, for the case in which both faces of the 
plate are equally cooled. In addition to exact formulas for the variable- 
conductivity plates, simpler approximate formulas are derived, and 
criteria for their applicability are indicated. 

2. It is. shown that the plate surfaces are always under the greatest 
tension, and the midplane is under the greatest compression, of all 
planes in the plate. 

3. A criterion is indicated for determining the surface cooling 
conditions under which the transient thermal shocks at the surface and 
midplane will be smaller than, equal to, or greater than the steady- 
state thermal stresses at those planes. The shocks exceed the steady- 
state stresses only when the initial surface-cooling conditions of the 

' transient state are more severe than the surface-cooling conditions of 
the steady state. The criterion is expressed for heat -generating plates 
by a simple relation between a steady-state heat-transfer parameter of 
the heat -generating plate and the Biot number of the transient state. 
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4. The dimensionless parameters governing the transient temperatures 
and. thermal stresses in materials of linearly varying conductivity are 
derived hy a similarity study of the conduction equation and boundary 
conditions of the transient state. A numerical technique for solving 
the transient-state equations is indicated in detail. The method is 
employed to obtain numerical values of the transient temperatures and 
stresses for a variety of parameter combinations, selected to test the 
thermal shock criterion. The results of the calculations substantiate 
the shock criterion. 

5. As a result of the study, an insight is obtained into the steady 
and transient temperature and stress mechanisms in heat-generating plates 
of constant conductivity and of conductivity that varies linearly -with 
temperature . 

Lewis Flight Propulsion Laboratory 

National Advisory Co mmi ttee for Aeronautics 
Cleveland, Ohio, April 21, 1953 
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APPENDIX A 
SYMBOLS 

The folio-wing symbols are used in this report: 

A symbol for the quantity 2 (K 8 /K s ^ Q ) (K 8 /K s ^ Q - l)(q'" L 2 /2K e 0 s,o) 

c specific heat of plate material, Btu/lb-°F 

E modulus of elasticity, lb/sq in. 

h heat transfer coefficient in transient state, Btu/sec-ft 2 -°F 

h st heat transfer coefficient in steady state, Btu/see-ft 2 -°F 

K thermal conductivity of plate material, Btu/sec-ft-°F 

K e thermal conductivity of plate material at temperature T e , 

Btu/sec-ft-°F 

K g thermal conductivity of plate material at surface temperature, 

Btu/sec-ft-°F 

L half -thickness of plate, ft 

m coefficient in equation for thermal conductivity, (K = n£T + n) 

N number of thin segments into which plate thickness is divided 

n constant in equation for thermal conductivity, (K = mff + n) 

q M heat flow per square foot of plate area normal to direction 

of heat flow, Btu/sec-ft 2 

q ,M heat generation per unit volume, Btu/sec-ft 3 

T temperature, °F 

x coordinate, distance of plane from plate midplane, ft 

a thermal diffusivity of plate material, ft 2 /sec 

(3 coefficient of linear thermal expansion, ft/ft-°F 

thickness of segment in a plate that has teen divided into 
N equal segments, ft 
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6 difference "between local plate temperature T, and constant 

reference temperature T e , 9 = T - T e , °F 

v Poisson's ratio 

p weight density of plate material, lb/ft^ 

a normal stress, lb/sq in. 

a* dime nsionless stress parameter, a* = cr(l-v)/|3E0 s ^ o 

x time, sec 

Parameters : 

Biot number 

Biot number with plate conductivity evaluated at temperature 
Fourier number 

Fourier number with plate conductivity evaluated at temperature 

T e 


hL/K 

hL/K e 

arr/li 2 

OeT/L 2 


Subscripts : 


av 

c 

e 

f 

i 


iil/2 

i 


m 


average 

coolant of the steady state 
envirpnment of the transient state 
fictive layer 

integral, multiple of unit distance (Ax), i = 1, 2, 3, . . ., 
and designates distance of local plane from prescribed' 
reference plane 

multiple of unit distance, (Ax) 

integral multiple of unit time -increment (At) , j = 1, 2, 3, 

midplane of plate 

maxi mum value taken on "by variable 


max 
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min- minimum value taken on "by variable 

o steady state 

s surface of plate 

y direction 

z direction 
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APPENDIX B 


TEMPERATURES 

The temperatures in a plate are governed by the conservation of 
heat in the plate interior and at the surface. Heat conservation in the 
interior of a heat -generating plate of variable conductivity is expressed 
hy the familiar equation 


pc 


(I) 


a ( k dT 

5 ^ ) 5x 


) 


+ q' 


(Bl) 


If T e is an arbitrarily assignable constant, and 9 is defined as 

(B2) 

Heat conservation at a surface where the plate transfers heat to a 
coolant is expressed by the standard relation 


9 = T - T e 

equation (Bl) may be rewritten 


1 


t I 



(B4) 


where 


0 s “ 


T - 
x s 


T c 


(B5) 


and 


0 c = T c - T e 


(B6) 


2920 



2920 


MCA TN 2988 


31 


In the following discussion the plate is assumed to he eq ually cooled on 
both faces, in both the steady and transient states. 


Steady State 

Heat -generating plate of constant conductivity. - In the steady state 
the left member of equation (B3) vanishes, and if the plate conductivity 
is constant, equation (B3) simplifies to 



-1 


1 1 T 


(B7) 


with the condition q"' > 0 because the plate generates heat. Direct 
integration of equation (B7) and successive application of the conditions 

that at x = 0, ^ =s 0, and at x ** ±L, 0 = Q G}0 > lead to the familiar 

form 



or 


6 


0 


s,o 


1 + 


q’” 1? 

2K0 s ,o 



(B9) 


The difference between local and surface temperatures is 



Heat-generating plates of variable conductivity. - If the plate 
conductivity varies linearly with temperature, the equation for conduc- 
tivity may be written 


Kt=mT + n= m0 + (mT e + n) = m0 + Kg 


(BIO) 
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Then equation (B3) "becomes for the steady state 


_d_ 

dx 


[(" + I] - 


(Bll) 


Direct integration yields 


d0 


(m0 + Kg) ^ = -q' ” x + Ci 


Inasmuch as = 0 at x == 0, Ci = 0. A second integration yields 
dx 


i me 2 + Kg0 = - 5-' — x2 


2 + C 2 


At the surfaces x = iL, 0 = 0 SjO > therefore. 


c 2 = ^ m0 s,o + Ke 0 s,o + — — g 


q ' " L 2 


Thus 



+ KgS -{±*£> 1,0 


+ KgSg^o 




0 (B12) 


From equation (BIO), 


Kg,o — m®s,o + Ke 

It is now assumed that Kg ^ Kg^ 0 , or ^s,o/^e / which is equivalent 
to assuming that the arbitrarily prescribable reference temperature T e 
is taken different from the steady-state surface temperature T S ^ Q , and 
the quantity 0 S ^ O = (T s ^ 0 - T e ) ^ 0. This assumption is made to permit 
solution for the ratio 0/0 g ^ o , and is not necessary if 0 g ^ o is not 
desired to enter in the denominator of the solution. Then 
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m ;= 


K, 


s, o 


K a 


8,0 


(B13) 


IO 

c!> 


O 


When equation (B13) is substituted into (B12) and solved for Q/8 s ,o ^ 
the quadratic formula, the result is 



Thus 0/0 s ^ Q in the steady state is governed by x/L, q'" L^/2K e 0 s ^ o , 
and K s ^ 0 /K e , and it is recalled that equation (B14) has meaning when 
q ,M ^ 0, and K g 0 /K e ^ 1. Equation (B14) is not inconvenient to use; 

nevertheless an approximate relation is instructive, and is discussed as 
follows : 

Approximate formula for 0/0 s in variable-conductivity plates. - 
An approximate formula for 8f®s,o noay he obtained from equation (B14) 
as follows. The equation is of the form 



On restricting considerations to values of $ > 0 (i.e., plates for 

which K decreases as temperature increases), on approximating \/ 1 - 4> 
by (l - <j>/2), and on substituting for $ in equation (B15), 
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(B17) 


On carrying through the operations in the right member of equation (B17), 


- 1 


s,o 




(B18) 


Alternately, 


0 


0 


8,0 


1 


+ 


q rtt L 2 
^s,o 0 s,o 



(B19) 


Equation (B18) shows that for the ranges of K s ^ 0 /K e and 4> in 
which equation (B18) is a satisfactory approximation, ( 0 /®s,o) “ 1 
varies directly with q ,M L 2 /2K e 0 s ^ o and inversely with K s ^ 0 /K e , and 
that the shape of 0/0 SjO with x/L is a parabola. Equation (B19), 

which is simply a reduced form of equation (B18), is of the s a m e form as 
the equation for a plate of constant conductivity, the conductivity 
being that at the temperature of the plate surface. Thus equation (B19) 
contains no provision for conductivity variation in the plate interior, 
and if the conductivity decreases as temperature increases, equation (B19) 
will underestimate the plate interior temperatures (while if the conduc- 
tivity increases with temperature the equation will overest im ate the 
interior temperatures ) . Equation (B19) also fails to indicate that a 
limit exists for the feasible value of q ,M L 2 /2Kg0 s Q if conductivity 

decreases as temperature increases, as is shown in the next section. 

In order to consider the range of 4> in which equation (B19) is a 
suitable approximation, equation (B15) may first be written 
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= 1 + 




^S.O 


K, 


s,o 


(B20) 


Ke 


A comparison of equations (B18) and (B20) shows that the approximation 
resides in the quantity r^--n , which is equal to ^ s ^o \ and 


r s _ - 

p 

l! . which is eaual to ^ a >o\ 

L 0 S,O 

L_ 

CD 

o 


measures the temperature rise in the plate interior. On replacing 
f T-T t 

1, on denoting the values in equations (B20) and 




■by 


-s.o 


e 


s.o 


(B17) by the respective subscripts "exact" and "approximate," and on 
dividing equation (B20) by (B17), the result is 


(®“^s,o) exact 
(T-T s }Q ) approximate 


= | (1 - \/l~ *) 

V 


(B21) 


The function of 4> in equation (B2l) has the limit unity as 4> becomes 
zero (x/L ->±l) , and increases steadily to the value 2 as <f> increases 
to its upper li m it (<j> = l) . Thus whereas equation (B19) necessarily 
gives the correct surface temperature, it underestimates the plate mid- 
plane temperature by the greatest amount, and the exact difference between 
midplane and surface temperatures may be twice as great as the value 
given by equation (B19) if $ is unity at the midplane. (From eq. (B16) 
it may be seen that 4> increases in the direction of the plate interior.) 
If the error in (T-T s ^ 0 ) is defined by 6 in the relation 


^s,o^ exact 

(T-T 8 ,o) 


approximate 


■d +£ ) = f 


V l - 4>) 


(B22) 


it is seen on solving for $ that the error s occurs when 


* 


46 

(1 + e )2 


(B23) 


THUS, for example, (T-T S)0 ) ej . act - 1.1 C^s, approximate vheD 

= 0.33, and (T— T a , 0 ) exact “ 1*^ (-C"-^s,o) approximate w ^en = 0.56. 
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Limitation on q MI L2/2Ke0 s ^ o when conductivity decreases with 
temperature. - Inasmuch as Q/9 s , q is always read, the expression under 

the radical in equation (B14) must never he less than zero. If the plate 
conductivity decreases as temperature increases a value of K s>0 /Ke 
smaller than unity is associated with a positive value of 0 g Q , and a 
value of Kg^ 0 /K e greater than unity is associated with a negative value 
of ®s,o* (The possibility of a negative 0 S ^ O arises because T e is 
arbitrarily assignable.) For either combination of Kg^ 0 /K e and 0 g Q , 

the quantity 2 vr— l) ~~ 1 is positive. On imposing 

v^s , o / \-“-s , o / \2K ePs,o) 

the requirement that the radicand never be smaller than zero, the follow- 
ing limitation on the feasible value of q'*' L 2 /2K e 0 g Q is obtained. 


q'" L 2 
^e\ 6 s,o\ 



(B24) 


where the symbol | | indicates absolute value. Physically the limitation 
arises because if the assumed law of conductivity variation is valid the 
conductivity becomes "zero" at the plate midplane when the equality in 
expression (B24) is satisfied, and becomes "negative" when q’ ' ' I^/SKe^s,© 

exceeds the right member of relation (B24). (if the plate conductivity 
increases as temperature increases there is no limitation on 
q«n x. 2 y r 2K e 0 s 0 ; the sole requirement is that the plate surface tempera- 
ture never acquire the value for which the conductivity is "zero.") 

Relation between q' 1 ' L^/2Ke0 s ,o and steady-state cooling condi- 
tions^ - The fact that a plate generating q’" uniformly throughout its 
interior is maintained at a surface temperature 0 SjO in the steady 

state implies definite cooling conditions in this state. This fact may 
be used to obtain a relation for q' *' L^'/2K e 0 g ^ o . From the heat balance 


*st ( e s,o - 9 C ) * 4"' L 
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Both sides of the equation may he m u ltiplied hy 


L/ 2Kee S)0 ; hence. 


g 1 " L 2 1 ( h st L \ ( 9 s,o ' 6 c\ 

2Ke^s,o Ke J\ ®s,o ) 


(B25) 
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APPENDIX C 

APPROXIMATE FORMULAS FOR STEADY-STATE STRESS IN HEAT- 
GENERATING PLATES OF LINEARLY VARYING CONDUCTIVITY 

The formulas derived in the Analysis section for surface and mid- 
plane stresses in heat -generating plates of linearly varying conductivity 
are 


and 


■where 


Ks, o 

l- 1 ‘ A ln 

( i +n/a\ 

Ke 

2n/A 

VI -n/A /_ 


2 ( i Ka ’°\ 



[ Ke ) 



* 


Ks,o 

Ke 

1 + 

1 - A ln( 

2\Ta 

f l+rJT\ 

V.1 

- ^ 

^1 - A 

2 

( , 

V *e j 




(Cl) 


(C2) 


(C3) 


Approximate formulas useful for rapid estimate of the stresses may he 
obtained as follows: From reference 4 (p. 9l), for A < 1, 


In f 1 + = 2 

>2 + + aV* + ... 

1 

H 

3 5 J 


(C4) 


If A 
beyond 


is substantially smaller than unity 


A 3 / 2 

■= may be neglected; hence. 


(e.g.. 


A < 0.2), terms 
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2a/T ^1 + 



Then 


2 VA ^(l 


(1 - A) 





On substituting equation (C6) into equation (Cl), 



(05) 


(C6) 


Replacing 2A/3 by its equivalent from equation (C3), and carrying 
through the algebraic operations, yields 



2 A'" LVaW>.,oY 1 + A\ 2 

3 ^ K 0;O /Ke J[ 2) 3 


f q"'L 2 \ 

[2Ks,O 0 8,O I 



(C7 ) 


The midplane stress may he approximated as follows: 
's/I - A by fl - and on using equation (C6), 


On approximating 



2 I 1 ’ f) - t - A > 


Thus 


* 

°m,o 


1 

3 



L 2 /2Ke e s,oA 
*s,o/*e / 


(1 - A) 


1 / q" 1 L 2 \ 

3 K,o fl 8,o J 


(1 


A) 


(C8) 
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It is seen from the first forms of approximations (C7) and (C8) that if 
A is very small by comparison ■with unity, the surface and midplane 
stresses are directly proportional to q ,M and inversely 

proportional to K^o/Ee* the second forms of equations (C7) and (C8) 
show that when A is very small by compar ison with unity, the plate 
behaves as one having a conductivity that is constant at the value of 
the plate surface. 

Comparison of equations (C7) and (C8) with formulas (13) and (14) of 
the Analysis for plates of constant conductivity shows that in a variable 
conductivity plate the surface temperature is farther from the average, 
and the midplane temperature closer to the average, than in plates of 
constant conductivity. 
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APPENDIX D 


CALCULATION PROCEDURE 

Temperature . - Steady-state temperatures were obtained by use of the 
equations derived in appendix B and given as equations (9) and (10) of 
the Analysis section. The transient temperatures were obtained as follows: 

The plate transient temperatures are determined by the initial tem- 
perature distribution and by the conservation of heat. Heat conservation 
in the plate interior and at the plate surfaces is expressed alternately 
by the well-known equations 

bd S ( v ~dd\ , * 

pc §? = V K (D1) 


and 

" Ks (IDs = h0s (D2) 


Inasmuch as only the early stages of the transient state were of interest 
a finite-difference solution, which for extended time spans might be less 
practical, was employed, as follows: In the accompanying sketch the 

designations i-1, i, i+1 refer to stations distant Ax from one another 
in the plate interior j i-l/2 and i+l/2 represent half-way stations. 

The curve ABC represents the local 9 distribution after jAT time has 
elapsed. If station i is considered in the next time interval At 


9 
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52 "1 *T TD 1 

0 i,j+i - e i,ji r 

e i,j+i - 0 i,,i] 

Wy 

J Llli j 

At-» 0 L 

AT J L 

At J 


(D3) 


If stations i-l/2 and i+l/2 are considered at time jAr the slopes 
of the 9 curve at these stations may he approximated hy the slopes of 
the chords AB and BCj thus. 






Ax 


. , fde \ 

i M Bx ) i 
i+ 2^ V° X /i+|,j 


(K 


i JL" 


e i+l,,3 - e i 3 




Ax 


(D4) 


(D5) 



the neglect of the limit process, and use of equations (D4) and (D5), 
leads to the approximation 


_a_ 

dx 




e i+i,j) - (*i-|„j + * 1 + 1 , 5 ) e !,3 + (*i-|,d ®i-i,d) 


(D6) 

X 


On substitution of equations (D3) and (D6) into equation (Dl), and on 
solution for 0^ j+1^ obtained relation is 



2920 



CG-6 'back 0262 


NACA TIT 2988 


43 


If the 0 dlstrihution at time jAr is known, equation (D7) permits 
calculation of the 0 distribution in the plate interior at time 
(j+l) At. 

The temperatures prescribed by equation (D7) are constrained by 
equation (D2), to which the fictive-layer procedure was applied as recom- 
mended in reference 5. In this procedure the assumption is made that 
after the initial instant the temperatures at the plate surface and at 
the interior station distant Ax/2 from the surface lie on a straight 
line having the slope indicated by equation (D2). It Is also assumed 
that this straight line may be extended a distance Ax/2 to the exterior 
of the plate, and that the so-introduced fictive layer of thickness Ax/2 
may be treated as a segment of the plate. If 9f is the 9 value at 

the edge of the fictive layer (fig. 10 ) the algebraic equivalent of the 
assumption is that 

, ( 9 t + 9 Ax/2\ 

(0f - 9/sx/2h f%9\ he s, j h \ 2 /j 

* = “ 


where 0^/ 2 is the value at station 1 of figure 10(a) or 



hAx 

0 f,j “ 0 Ax/2,J 


i+ ““ 


2Ks,j_ 


(D8) 


Whereas this relation is employed for 0f after the initial instant, 
the value of 0^ at zero time is assumed to be given by the relation 


( 0 f - 0 s)o 
Ax/2 



h0 


8 ,0 


Ks,< 


or 


0 f,o “ 0 s,o 



(D9) 
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A knowledge of Qf permits use of equation (D7) to determine temperatures 
in the plate interior. 

In applying equations (D7) through (D9) the increments At and Ax 
were constant throughout the calculation at values determined hy the con- 
vergence criterion indicated in reference 6. 


At 

( A *) 2 





(DIO) 


Inasmuch as Ac and Ax may he chosen arbitrarily provided they satisfy 
equation (DIO), the choice in this analysis was made according to the 
relation 


A'E 1 

(Ax)^ /•^max'\ / hAx A 

V P c ) \ ^min ) 

wherein K max and n were the respective maximum and minimum values 

of plate conductivity occurring during the total time span under consid- 
eration for thermal shock. This choice insured that equation (DIO) would 
always be satisfied. 

Dimensionless parameters for presentation of specific numerical 
results in a more generally useful form were derived by Nusselt 1 s method 
of applying the similarity principle to heat transfer, as described in 
reference 7. The derivation of heat flow parameters is shown in 
appendix E. 

Stress. - The stress was computed by use of equation (8) of the 
Analysis section. In the steady state, equations (13) and (14), and 
(17) and (18) of the Analysis were employed. In the transient state, 
the quantity (0/0 s ^ o ) av was obtained by planimeter integration under 

the transient curve of (®/® s 0 ) over the plate thickness. 
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APPENDIX E 

PARAMETERS THAT GOVERN TRANSIENT TEMPERATURE IN PLATES 
OF LINEAR LI VARYING CONDUCTIVITY 

For the transient state it is assumed that the plate suddenly ceases 
to generate heat and is simultaneously Immersed in a cooling environment 
of constant temperature T e , to which the plate transfers heat hy con- 
stant heat-transfer coefficient h. Equation (B3) becomes 

pc * JlYk (El) 

d'C bx V bx J 

and equation (B4) becomes 

- he s (E2) 


It may he noted that whereas in the steady state the temperature 
T e was arbitrary, the use of 0 in equation (E2) implies that for the 

transient state T e must be equal to the temperature of the transient- 
state cooling environment. In practice it is convenient to use the same 
value of Te for both the steady and transient states, because the 
results of the transient-state calculations depend on a knowledge of the 
steady-state temperature distribution, and it is convenient to have the 
original temperature distribution in the same terms as employed for the 
transient-state calculations. If the thermal conductivity of the plate 
is sufficiently linear in the temperature range of interest and can be 
written 


E = mT + n 


(E3) 


it can be rewritten as 


K = m(T - T e ) + (mT e + n) = + ^ 


(E3a) 
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Equations 


(El) and (E2) "become in terms of 


0 


and 



(Ela) 



(E2a) 


The dimensionless groups which govern the transient temperatures 
may he determined hy Nusselt 's method of applying the similarity prin- 
ciple to heat transfer (ref. 7). Accordingly, plate systems 1 and 2 
are now considered, and are required to he similar in geometry, time, 
and temperature. If x, t, and 0 represent local values, and if a, 
h, .and f are constant independent of x, t, and 0, the similarity 
requirement is expressible as 

h. 

x*j^ ~ axg ~ Lg x 2 


T! = trcr 2 (E4) 


0]_ =S f0g 


The respective system constants, pc. Kg, h, and m hear the constant 
ratios 
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(pc)l = g(pc) 2 
K e,l = i Ke,2 
h l = i *2 


Hi = l 


mg 


(E4a) 


It is noteworthy that m i 4 “2 implies nonparallelism of K with 0 

in systems 1 and 2. The temperatures in plates 1 and 2 satisfy equa- 
tion (Ela), i.e.j 


(p?)i 





*e,l 



(pc) 2 


^2 ( Q ^ 2 \ K h2 °2 

^2 Sx 2 \ 2 Sx 2/ ,2 ^x 2 2 


(E5) 


(E6) 


When equations (E4) 


and (E4a) are employed in equation (E5) 



Equations (E5a) and (E6) are identical if 



(E7) 


(#)©•* 


(E7a) 
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Use of equations (E4) and (E4a) in equation (E7) gives 


ga 2 _ (P c )l x l 2 ^2 K e,2 
~ (pc)2 x 2 2 t-l Ke,l 



(E8) 


Thus, one dimensionless group is Fourier's modulus "with the plate thermal 
diffusivity evaluated at the environment temperature. Equation (E8) 
shows that the times at which the temperatures of similarly located points 


sire to he compared are in the ratio 
and (E4a) are used in equation (E7a) 



^1 


When equations (E4) 


If _ m l 9 l ^e,2 
1 = “2 0 2 *e,l 


or 


fl = *^,1 f2 

e 2 ^,2 m l 


(E9) 


Equation (E9) shows the actual ratio existing between temperature dif- 
ferences 0-^ and 0g of similarly located points when considered at 

the corresponding times designated by equation (E8). The explicit value 
of the 0 ratio does not appear in the case of materials of constant 
thermal conductivity. 

The dimensionless group that must have the same value in the two 
plate systems if equation (E9) is to be satisfied may be obtained by 
considering the surfaces of the two plates at zero time. If 0 g Q rep- 

sents initial t emp erature difference between surface and environment , 
equation (9) requires that at t =* 0, 
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But from equation (E3a), m0 s ^ o = Kg^ 0 “ where K s 0 is the plate 
conductivity at the initial surface temperature. . Hence, 


t- 



or 




(E10) 


Thus for two plates of linearly varying conductivity, temperature similar- 
ity can exist if the ratio of plate conductivity at the initial surface 
temperature to the plate conductivity at the environment temperature 
is the same in hoth systems. 

For systems in which equation (E10) is satisfied equation (E4) is 
also valid and 


fl - K e,l ^2 ( e s,o)l 

°2 ^e,2 “l ( 9 b , o ^2 


and 


fc)i _ fe)a 


(Ell) 


Systems 1 and 2 hoth satisfy equation (E2a) . Thus 


(h0 s )l = 



(E12) 
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and 


(M s ) 2 - 



fbe\ 


r f%o\~ 

— — 

h yj 

2 



(E13) 


On using equations (E4) and (E4a) in equation (E12), 


Jf(Ms)2 




(E12a) 


Equations (E12a) and (E13) are identical if 


(jf) (a/if) «. 1 


and 



On using equations (E4) and (E4a) in (E14), 


ja _ h l h. K e,2 
1 *2 h K e,l 


(E14) 


(E15) 


or 



(E16) 


Biot's modulus with, plate conductivity evaluated at environment tempera- 
ture is another governing parameter. For temperature similarity of two 
plate systems the heat-transfer coefficients must he in the ratio 
(L/K e ) 2 /(L/K e )i. 
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Equation (E15) yields again the requirement (ia0/K e )^ = (mS/K^g, 
which was obtained from equation (E7a). 

Equations (E4), (8), (10), (ll), and (16) show that results of 
numerical calculations for any type of initial temperature distribution 
can be presented in a form useful for similar systems having the same 
type of initial distribution by plots of 9/Q q , 0 against hl/Kg, with 

a e T/L 2 as parameter, for prescribed values of x/L and Kg ^ 0 /E e . For 

general usefulness, dimensionless temperatures for the entire range of 
probable K s ^ 0 /K e values, and for various values of hl/Kg and pres- 
cribed initial distribution, require calculation, (numerical results 
for Kg 0 /K e = 1 and 0.15 are presented in this report for the initial 

temperature distribution occurring in a plate which in the steady state 
generates heat uniformly throughout its interior and is equally cooled 
on both faces . ) 

It is noteworthy that the dimensionless groups for plates of linear 
conductivity and constant volumetric specific heat, exchanging heat in 
the transient state with an environment of constant temperature T e , 
all require the plate conductivity to be evaluated at the environment 
or surface temperature and include the dimensionless group K s 0 /K e , 

the ratio of plate conductivities at initial surface and at environment 
temperatures . The ratio K g } 0 /K e measures the rapidity with which the 

plate conductivity varies with temperature when 9 8j0 and K e are 
known, or indicates the AK experienced in going from T e to T g . 
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Figure 1. - System under consideration in steady state. 
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Figure 2. - Distribution, of for several values of Eg 0 /B^ and q.' ' , L 2 /2K e 0 0 Q . 
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Dimensionless temperature parameter, (q 1 ' ' L^/^KgSg^) 


Figure 3. - DimenBlonless surface and midplane stresses In heat-generating plates of constant 
conductivity and of conductivity that varies linearly with temperature, as functions of 
dimensionless temperature parameter q ' 1 'L^/aKgSg 0 _ and conductivity parameter Kg o/Kg. 
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0 .2 .4 .6 .8 1.0 

Location, x/L 


(a) Temperature distribution. 



0 .1 .2 .3 .4 .5 

Fourier number, ax Ju 


(b) Surface and midplane stresses. 

Figure 5. - Transient temperatures and stresses in formerly-beat-generating 

plates of constant conductivity, q. ' ' 'L 2 /2K0 S ^ O = 0.01, KL/Kq = 0 . 01 < 2 • 
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0 .2 .4 .6 .8 1.0 

Location, x/L 


(a) Temperature distribution. 



(b) Surface and midplane stresses. 

Figure 6. - Transient temperatures and stresses in formerly-heat-generating 

plates of constant conductivity, q.' ' 'L 2 /2K0 S ^ O <= 0.01, hL/t = 0 . 03 > 2 • 
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Figure 7 


(t) Surface end midplane stresses, 03 and ctJ. 

Transient temperatures and stresses In formerly-heat-generatlng plates of 


yuj. CJ I . _ A1CUIOJ.UUU w wiUf' VO. ^ — r " — — - 

variable conductivity. ^ b } q/^o “ l'" L ^/ K e 0 s,o = 0 . 01 , JiL/Eg ° ° °* 02 ‘ 




60 


NACA TBT 2988 



0 .2 .4 .6 .8 1.0 

Location, i/l 


(a) Temperature distribution. 



Figure 8. - Transient temperatures and stresses in f ormerly-heat -generating plates of variable 
conductivity. K g>0 /K s - 0.151, q" , L 2 /2E 0 e BjO - 0.01, hL/K, = 0.01<2 " 


Jrig plates 
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(a) Temperature distribution. 
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(a) Plate divisions and fictive layer. 



(t) Assumed fictive -layer temperature 6 S ^ 0 at start of transient 



(c) Assumed fictive layer temperature at time JAx of 


transient state 


- ^ m f - c]/! 1 + Cj} 


L L “'-J 

Figure 10. - Details of finite -difference procedure for obtaining temperatures 
in transient state. 
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